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Abstract. We report in this paper the numerical simulations of the capture into the 3:1 mean- 
motion resonance between the planet b and c in the 55 Cancri system. The results show that 
this resonance can be obtained by a differential planetary migration. The moderate initial eccen- 
tricities, relatively slower migration and suitable eccentricity damping rate increase significantly 
the probability of being trapped in this resonance. Otherwise, the system crosses the 3:1 com- 
mensurability avoiding resonance capture, to be eventually captured into a 2:1 resonance or 
some other higher-order resonances. After the resonance capture, the system could jump from 
one orbital configuration to another one if the migration continues, making a large region of 
the configuration space accessible for a resonance system. These investigations help us under- 
stand the diversity of resonance configurations and put some constrains on the early dynamical 
evolution of orbits in the extra-solar planetary systems. 
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1. Introduction 

Up to date, more than 200 extra-solar planetary systems have been found. Those host- 
ing more than one planet are multiple planet systems. In some of the multiple planet 
systems, planets are observed to be locked in mean-motion resonance (MMR), for exam- 
ple, the well-known 2:1 MMR in GJ876 system (e.g. Marcy et al. 2001) and the 3:1 MMR 
in 55 Cancri system (e.g. Zhou et al. 2004). In a 3:1 MMR, at least one of the three reso- 
nant angles (6*1 = Ai — 3A2 -l-2ci7i, 6*2 — Ai — 3A2 -I- 2n72, 6*3 — Ai — 3A2 +072) librates, 
where Ai^2 and wi^2 are the mean longitudes and periastron longitudes of the inner and 
outer planet respectively. The resonant systems are particularly attractive because not 
only of the complicated dynamics of the resonance but also of interesting information 
about its origin and early evolution buried in the configuration and dynamics. 

The 55 Cancri system is the only example of 3:1 MMR found in the extra-solar plan- 
etary systems till now. Four planets have been reported in this system and two of them, 
55 Cnc b and 55 Cnc c, seem to be locked in a 3:1 MMR. In our previous work, we have 
found that they are most likely in one of the three possible configurations (Zhou et al. 
2004). G. Marcy declared in his lecture in this Symposium that a fifth planet has been 
dug out from observing data of this system, enriching the whole story with more conno- 
tations. In the dynamical simulations of their 5-planet solution (Fischer et al. 2007), they 
did not found the resonance between the two planets, although their orbital periods are 
very close to the 3:1 commensurability. However, the current orbit determinations may 
not be robust enough to determine the resonant angles, especially for orbits with low 
eccentricities as be given in this new determination, thus the resonant character of this 
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Table 1. Orbital elements and masses. The mass of the central star is l.OSM© (McArthur et 
al. 2004). In this paper we adopt the planet masses by assuming sini = 1. 



Parameter 


Caneri c 


Caneri b 


Caneri e 


Caneri d 


Msini (A/j) 


0.045 


0.784 


0.217 


3.92 


P (days) 


2.81 


14.67 


43.93 


4517.4 


a (AU) 


0.038 


0.115 


0.240 


5.257 


e 


0.174 


0.0197 


0.44 


0.327 



pair of orbits may change as more observations are included (Beauge 2008). Furthermore, 
what we will discuss in this paper has some general significance so that we may still adopt 
the orbital elements and planetary masses from the previous literature (McArthur et al. 
2004), as listed in Tabled 

It is widely accepted that a planet may experience (generally inward) migration, and 
through differential migrations of planets, different commensurabilities among planets 
can be attained (e.g. Nelson & Papaloizou 2002). For example, the 2:1 MMR in GJ876, 
whose configuration can be obtained through the inward migration of the outer planet 
(Lee & Peak 2002). 

On the other hand, the resonant capture probability depends on the order of the 
resonance, the migration rate and the initial planetary eccentricity (Quillen 2006). The 
formation process of the 3:1 MMR may be very different from the one of the 2:1 MMR. 
As the only example of 3:1 MMR in extra-solar planetary systems, the 55 Caneri system 
deserves a detailed investigation. In this paper, we will report our numerical simulations 
of the capture into resonance, and the evolution of the resonance thereafter, provided 
the migration continues. 

2. Numerical simulations of resonance formation 

There are more planets than the two (planet b and c) in our focus of attention in this 
planetary system. However, besides these two planets, the massive one (55 Cnc d) is too 
far and the close one (55 Cnc e) is too small. Moreover, even the newly reported planet 
(55 Cnc f) is neither very massive (M/sini — 0.144Mjup) nor very close to this planet 
pair (a/ = 0.781 AU, P/ = 260 days). Therefore as a simplified model, it is reasonable to 
discuss only the two planets under the influences of the central star and the disc. 

2.1. Numerical Model 

As usual, the influence of the disc on the planets is simply simulated by an artiflcial force, 
which is acting on the outer planet to drive it to migrate inward. The force is defined so 
that the semi-major axis of the planetary orbit will change following an exponential law: 

a(i) ao + Aa X e"*/^. (2.1) 

where ag is the initial semi-major axis of the planet, and t is the timescale of the mi- 
gration. The two planets are assumed to be located initially on orbits with semi-major 
axes of 1.25 AU and 0.5 AU, with an initial semi-major axes ratio of 2.5 (periods ratio 
3.95), which will evolve downward to the value corresponding to the 3:1 MMR (2.08). 
Obviously the value of t determines the migration speed. The smaller the r, the faster 
the migration, and vice versa. Different r values were adopted in different papers. The 
standard migration scenario (type II) gives an estimate of migration rate (Ward 1997): 
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Taking the typical parameters a and H/a, the r for a planet starting from 1.25 AU is 
^ 1.43 X lO^yrs. Some other rules have been applied in different papers. According to 
these rules, the r would be, for example, 1.56 x lO^yrs (Lcc & Peale 2002), 1.38 x lO^yrs 
(Nelson & Papaloizou 2001), and 2.65 x lO^^yrs (Kley 2003; Kley et al. 2004). In our 
simulations, it is set to be 2 x 10^,1 x 10^,2 x 10^ and 5 x lO^yrs. Some large values 
(meaning slow migration) arc adopted here to guarantee the capture into the 3:1 MMR. 

Although it is still not very clear now how a certain disk will damp or stimulate the 
eccentricity of an embedded planet, people hope the disk will circularize an orbit if it was 
too eccentric. In order to control the unlikely eccentricity increasing during the artificially 
induced migration, we have included an eccentricity-damping force in our simulations. A 
parameter K was introduced to describe the damping rate, with which we control the 
eccentricity e of a planet by: 

\e/e\= K\a/a\. (2.3) 

In our simulations, the parameter K has the values (no damping), 1, 10 and 100. 
The damping is always put on planet 55 Cnc c with the migrating force modifying its 
semimajor axis simultaneously. 

In the simulations, the two orbits are nearly coplanar (in the practice of numerical sim- 
ulations we adopt an arbitrarily defined small inclination between these two orbits, say 
0.1 degree), and the initial orbital eccentricities of planet 55 Cnc b (e^) and 55 Cnc c (e^) 
are set to be 0.001 (a typical value for a nearly circular orbit) as well as 0.01, 0.05, 0.1, 0.2. 
For each mode (with certain r, ii", e?, 6°), we numerically integrate 100 systems with ran- 
domly selected orbital angles (longitudes of pcriastrons, mean longitudes, and ascending 
nodes) . Each testing try is integrated up to a time of 2t. During an integration, if the two 
semi-major axes are locked in a definite value and keep this value at least for ten percent 
of the integration time, that is t/5, we say the system is in the given commensurability. 

2.2. Results 

Generally, the migration may drive the two planets into a given commensurability. The 
inner planet 55 Cnc b will also migrate after it has been captured by a commensura- 
bility. The final configuration of a system depends sensitively on the migration speed 
and the initial conditions. We summarize some remarkable outcomes from our numerical 
simulations as follow. 

(1) The migration with r ~ 2 x 10^ yrs is too fast to form the 3:1 MMR. If we start from 
nearly circular orbits (ej = 6° = 0.001) and neglect the eccentricity damping {K = 0), 
all the 200 runs for r = 2 x 10"' yrs and r = 1 x 10^ yrs cross the 3:1 MMR without being 
trapped and are eventually captured into the 2:1 resonance. The 3:1 MMR is observed 
only when t = 2 x 10^ yrs. When r = 2 x 10'* yrs, the number of simulation runs, in 
which the 3:1 MMR forms, is very small, if not null, no matter what ej, 62 and K are. 

(2) A slow migration favors the formation of the 3:1 MMR. For initially near circular 
orbits, when r = 2 x 10^ yrs, we have obtained 53 out of the 100 runs trapped in the 3:1 
MMR and the number increases to 100 when r = 5 x 10^ yrs. 

(3) The initial eccentricities strongly affects the resonance trapping. A nonzero but 
small initial eccentricity of cither orbits (ej or ~ 0.01) will increase significantly the 
probability of the 3:1 MMR. For the runs with e° or 62 = 0.01, although there is no 3:1 
MMR when r = 2 x 10^ yrs, we found 163 out of the totally 200 runs trapped in this 
resonance when r = 1 x 10^ yrs. Higher initial eccentricities (e^, ^ 0.05, 0.10) will bring 
in more higher-order resonances, such as the 5:2, 7:3 and 7:2 MMRs. It is interesting 
to see that a few runs were trapped in the 11:4 resonance. Even higher eccentricities 
(e?,e2 ^ 0.20) generally leads the system to catastrophic planetary scattering. 
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(4) The eccentricity damping affects the resonance capture too. Without the eccen- 
tricity damping {K — 0), higher initial eccentricities would likely lead to high order 
commensurabilities, while low initial eccentricities would cause eccentricities to grow too 
much after the planets are locked into a low-order resonance. On the other hand, a high 
K value restrains the eccentricity increasing and therefore causes the system to evolve 
preferably to a 2:1 MMR rather than the 3:1 MMR. For example, when K = 100 all the 
test systems starting from circular orbits (ej = 6° = 0.001) are driven into the 2:1 MMR 
if r ^ 2 X 10^ yrs. Even in a very slow migration (t = 5 x 10^ yrs), the 2:1 MMR captures 
41 out of the 100 runs. So we need a suitable damping rate for the formation of 3:1 MMR. 
In all our simulations, the most probable modes for the formation of the 3:1 MMR were 
a slow migration (t ^ 1 x 10^ yrs), moderate initial eccentricities (6^2 = 0.01 ^ 0.05) 
and moderate eccentricity damping {K 10). 

3. Evolution of orbital configuration in resonance 

After being captured into the 3:1 MMR, the system will continue to evolve if the 
migration does not halt. In this section we will discuss the evolution of the system in the 
resonance. 




Figure 1. Periodic solutions in the 3:1 MMR. (a) Periodic solution families on the (ei, 62) plane. 
The stable symmetrical ACR solutions are thick solid curves, the unstable symmetrical ACR 
solutions are short-dashed curves, and the thin solid curves indicate the asymmetrical solutions. 
The value of {61, Aro) of the symmetrical solutions appear as labels along the curves. We also 
label different solution families with A, B, C and D. The eccentricities (ei, 62) = (0.0197, 0.44) 
in Table [1] are indicated by a cross, (b) The variation of the angles 81 (dotted), 62 (dashed) and 
Auj (solid) along the solution family A. Note Aw begins from 180° (symmetrical ACR). (c) 
The same as (b) but for solution family C, where Azu reaches at symmetrical value 0° at the 
end. 

3.1. Periodic solutions in 3:1 MMR 

Fig.[T]shows the periodic solutions, both stable and unstable, symmetrical and asymmet- 
ical, computed using the mass ratio of the planets 55 Cnc b and 55 Cnc c and assuming 
coplanar orbits. (For other mass ratios and for details of the method used see Michtchenko 
et al. (2006)). 
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In the real system, if we assume that it evolved from nearly circular orbits through an 
adiabatical migration, as suggested in many references, the configuration of two planets 
captured into a 3:1 resonance should be located on or near the family A of periodic 
solutions shown in Fig.[T]Ja). In such a migration, the angles {0i,Ava) evolve from the 
symmetrical (tt, tt) values to asymmetrical values. Since each solution family is not con- 
nected by stable periodic solutions to other familie^, there is no reason to expect a 
system in resonance with orbital configuration far away from the family A. In fact, the 
eccentricities (ei, 62) in Table[T]are not on the family A. However, we have found that this 
configuration could be obtained by a resonance captured from initial eccentric orbits. An 
example of such kind of resonance trapping is shown in Fig.[2lja). Surely a perturbation 
happening to a system in the equilibrium periodic solution may also lead to this observed 
configuration. We just show here another possibility of the resonance formation. 

3.2. Jump between solution families 

If the disc does not disappear just after the resonance capture happens, the outer planet 
will continue to migrate inward accompanied by the inner planet locked in the resonance. 
During such after-capture migration stage, the system will evolve along the solution 
family A as shown in Fig.[l] provided the migration is "adiabatic" . But if the migration 
is faster (not adiabatic), as some of our numerical simulations, the solution may jump 
from a family to another. An example of jump from family A to C is illustrated in 
Fig. nub). Comparing the variations of angles in Fig.O^c) and Fig.[Hb), (c), it is clear 
that the jump happens around t = 1.9 x lO'^yrs. At this moment, if the system goes 
along family A, the angles di and Aw will go up and cross the value 180° one after 
another with the increasing eccentricity 62- But they fail and turn around following 
another evolving route, that is, family C. This turnaround is due to some perturbation 
(in our case it's the migration itself), and it happens around the critical point where the 
apsidal difference Atn crosses the symmetrical value 180° in an asymmetical ACR. In 
fact, the stable region of motion around this critical turning point (on family A with 
(ei, 62) = (0.0423, 0.623) in Fig.flja) ), is very small. Therefore a slight disturbance could 
drive the system out of the family A. 



(a) (b) (c) 




0.0 0,1 0.2 03 04 05 0.6 07 0,8 0.9 10 00 0.1 02 03 0.4 5 3 7 0.3 9 1.0 10000 20000 30000 40000 



e, Time (yrs) 

Figure 2. The evolution of orbital configuration in the 3:1 MMR. (a) An example of evolution 
starting from eccentric initial orbits. The eccentricities in Table[l]can be reached by the solution 
in this case. Solution families are shown in the same way as in Fig.[T] and the dots indicate the 
evolution of one run in our numerical simulations, (b) An example of orbital configuration jump 
from solution family A to family C. (c) The time variations of angles 6i (small dots), 62 (open 
circles) and Auj (solid thick curve) for the numerical simulation shown in (b). 

f A and B are not connected because they have different 9i, Azu values at the two crossings 
in Fig.ETa) 
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The jump from family A to family B has also been observed in our simulations. The 
angles Oi = 62 = 0° in family B, that is to say, both wi and vj2 have to adjust to 
new values in a jump from A to B. On the other hand, as we see from Fig.[T]^b),(c) and 
Fig-HUc), we need only an adjustment of Wi but not VJ2 for a jump from family A to 
C, since this jump does not affect the 62- When the jump happens, the inner planet is 
on a near circular orbit (ei = 0.0423) while the outer one is on a highly eccentric orbit 
(e2 = 0.623). As a result, it is much easier to adjust wi than vd2. This analysis tells why 
we observe much more A to C jumps than A to B jumps. 

These jumps between different solution families make it possible for a resonant system 
to occupy much larger potential volume in the orbital elements space, increasing the 
diversity of resonant configuration in extra-solar planetary systems. 

4. Conclusions 

The 3:1 mean-motion resonance in the 55 Cancri planetary system is far from a nec- 
essary result of the differential planetary migration. Our numerical simulations showed 
that the favourable scenario for the formation of the 3:1 resonance is moderate initial 
eccentricities (0.01 ~ 0.05), relatively slower migration (t ^ lO^yrs) and suitable ec- 
centricity damping rate {K ~ 10). After being captured in the resonance, the system 
may exhibit some evolutions different from the behaviours in an adiabatic migration, 
and evolve to some unexpectable orbital configurations. All these put some constraints 
on the planetary migration and early dynamical evolution in the extra-solar planetary 
systems. 
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